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ABSTRACT 

NaStl (aka Wolf-Rayet 122) is a peculiar emission-line star embedded in an extended 
nebula of [N ii] emission with a compact dusty core. This object was characterized by 
Crowther & Smith (1999) as a Wolf-Rayet (WR) star cloaked in an opaque nebula 
of CNO-processed material, perhaps analogous to rj Car and its Homunculus nebula, 
albeit with a hotter central source. To discern the morphology of the [N ii] nebula 
we performed narrowband imaging using the Hubble Space Telescope and Wide-held 
Camera 3. The images reveal that the nebula has a disk-like geometry tilted ssl2° 
from edge-on, composed of a bright central ellipsoid surrounded by a larger clumpy 
ring. Ground-based spectroscopy reveals radial velocity structure (±10kms“^) near 
the outer portions of the nebula’s major axis, which is likely to be the imprint of 
outhowing gas. Near-infrared adaptive-optics imaging with Magellan AO has resolved 
a compact ellipsoid of ATg-band emission aligned with the larger [N ii] nebula, which 
we suspect is the result of scattered Hei line emission (A2.06/im). Observations with 
the Chandra X-ray Observatory have revealed an X-ray point source at the core of 
the nebula that is heavily absorbed at energies < 1 keV and has properties consistent 
with WR stars and colliding-wind binaries. We suggest that NaStl is a WR binary 
embedded in an equatorial outflow that formed as the result of non-conservative mass 
transfer. NaStl thus appears to be a rare and important example of a stripped-envelope 
WR forming through binary interaction, caught in the brief Roche-lobe overflow phase. 

Key words: stars: individual (NaStl) - stars: emission-line - stars: peculiar - stars: 
Wolf-Rayet - stars: mass loss 


1 INTRODUCTION 

The physical processes by which massive stars shed their 
hydrogen envelopes and become Wolf-Rayet (WR) stars 
present an interesting challenge to the standard stellar evo¬ 
lutionary paradigm (see review by Smith 2014 and refer¬ 
ences therein). Mass loss via steady stellar winds is punc¬ 
tuated by discrete, episodic, and sometimes explosive mass 
ejections that occur during the post main-sequence, in the 
so-called Luminous Blue Variable (LBV) phase (Smith & 
Owocki 2006). The physics that trigger and drive these in¬ 
tense bouts of eruptive mass loss are still not understood, 
yet this mode is suspected to be a key channel by which 
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H-depleted WRs form. In addition to discrete eruptions, bi¬ 
nary mass transfer via Roche lobe overflow (RLOF) is an¬ 
other very important (perhaps even dominant) channel for 
producing stripped-envelope WRs, and may even be critical 
to account for Type Ib/c supernova progenitors (Smith et 
al. 2011a; Dessart et al. 2012; Hachinger et al. 2012). Since 
the time-scale of eruptive mass loss in such systems is brief, 
only a precious few Galactic examples exist, especially for 
cases where stellar multiplicity (i.e., binary vs. single) is an 
influential factor. A particularly interesting example is RY 
Scuti - a massive binary undergoing mass transfer that re¬ 
cently resulted in an intense bout of equatorial mass loss, 
and where one of the stellar components is likely on its way 
to becoming a stripped-envelope WR (Smith et al. 2011b). 
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NaStl, aka WR122, is an evolved massive star that 
has defied characterization. It was discovered by Nassau & 
Stephenson (1963), who proposed a WR classification based 
on the strong emission-line spectrum. It was later reclassified 
as a late nitrogen-type WR (WN type) based on nitrogen 
emission-line ratios in the optical spectrum (Massey & Conti 
1983), then later proposed to be a B[e], 0[e], or Ofpe/WN9 
star based on its optical-IR spectral properties (van der 
Hucht 1997). Subsequently, more detailed multi-wavelength 
investigation by Crowther & Smith (1999; hereafter CS99) 
showed that although NaStl is indeed a hot luminous ob¬ 
ject with L ~ 10®'^-10®’^ 1/0, its emission-line spectrum is 
mostly of nebular origin, not stellar. Moreover, ground-based 
narrowband imaging through an [N ll] filter centered near 
wavelength A6584 A revealed the presence of a compact cir- 
cumstellar nebula with a diameter of «7”, whose elongated 
yet poorly resolved dimensions gave the impression of pos¬ 
sible bipolar morphology. No comparable extended compo¬ 
nent was inferred for H l gas based on the much weaker de¬ 
tection of nebulosity in images through a broader Ha-|-[N ll] 
filter (which could simply be attributable to [N ll] emission 
in the passband). He l and He ll narrowband images revealed 
no extended component either. 


Derivation of the chemical abundances in NaStl’s neb¬ 
ula via quantitative spectroscopy by CS99 showed it to be 
extremely nitrogen enriched yet carbon and oxygen defi¬ 
cient, which indicates that the nebula consists of highly 
CNO-processed material ejected from an evolved massive 
star. CS99 pointed out that the only other massive star 
known to be enshrouded in a nebula having a compara¬ 
ble level of nitrogen enrichment is the LBV-I-OB binary rj 
Car and its Homunculus nebula (Davidson et al. 1986; Du- 
four et al. 1997; Smith & Morse 2004). However, CS99 note 
that the chemical abundance ratios in the NaStl nebula are 
characteristic of a WR, and probably inconsistent with a 
red supergiant or H-rich LBV. Moreover, the derived ion¬ 
ization potential requires that the central ionizing source 
be T » 30 kK, which suggests a spectral subtype earlier 
than WN6-7. The absence of the characteristic broad emis¬ 
sion lines from a WR wind, however, implies that the stel¬ 
lar system is deeply embedded in an opaque nebula. High- 
resolution spectroscopy by CS99 revealed that the nebular 
lines have a double-peaked morphology, which suggests com¬ 
plex geometry for the outflow. 

NaStl also exhibits strong thermal excess from circum- 
stellar dust, but from a much more compact region than the 
[N ll] nebula. CS99 presented evidence for a two-temperature 
spectral energy distribution in the infrared, consisting of 
warm (T ~ 700 K) and hot (T « 2000 K) components. 
A mid-infrared spectrum by Smith & Houck (2001) showed 
it to be mostly featureless in the 8-14 /rm wavelength range, 
with the exception of two narrow nebular lines of [Fe ll] and 
[Nell]. The results of mid-infrared interferometric imaging 
at 10 fim wavelength by Rajagopal et al. (2007) were con¬ 
sistent with at least two infrared components having angular 
sizes within 50-60 mas. Based on the likely distance range of 
1-3 kpc and on the measured outflow velocities from CS99, 
the results of Rajagopal et al. (2007) imply that the dust 


originates from within <100 AU of the central sourc^ and 
thus has a short dynamical timescale of only 30 yr, which 
implies only very recent and likely ongoing dust formation. 
By comparison, the outer [N ll] nebula has a much larger 
size of «7” (6,800-22,400 AU). The presence of an addi¬ 
tional component of faint thermal dust emission associated 
with the extended nebular structure, to our knowledge, has 
never been assessed with sufficient spatial resolution. 

The collective properties of NaStl led CS99 to con¬ 
clude that this object could be a WR cloaked in an rj Car¬ 
like ejection nebula, while the interferometry of Rajagopal 
et al. (2007) motivates a search for extended structure in 
the near-infrared. Here, we investigate the situation further 
by obtaining high-resolution narrowband images of NaStl’s 
nebula in the light of [N ii] using the Hubble Space Tele¬ 
scope (HST), and adaptive-optics (AO) imaging in the near- 
infrared with Magellan AO. Moreover, the possible presence 
of a WR star in an p Car-like system also warrants the 
search for an X-ray counterpart to NaStl, which could be 
expected to arise from shocks produced by hypersonic wind 
flows that are often associated with WRs. NaStl was not 
detected by the ROSAT all-sky Survey, but this is not par¬ 
ticularly surprising, given the relatively shallow sensitivity 
of this survey, in conjunction with the estimated distance 
to the source and the anticipated absorption to X-rays with 
energies < 1-2 keV from interstellar material. Thus, we have 
obtained deeper X-ray imaging observations of NaStl using 
the Chandra X-ray Observatory. The resulting data from 
these space-based and AO facilities provide new important 
clues for the nature and evolutionary status of this object. 


2 OBSERVATIONS AND RESULTS 

2.1 Narrowband [N ll] imaging with RST/WFCS 

NaStl was first observed with HST on 2013 March 20 
(UT 02:34) using the Wide-field and Planetary Camera 3 
(WFC3). Five 176s exposures were obtained through both 
the F658N and F645N filters, which have peak transmis¬ 
sion wavelengths centered on the [N ll] A6584 A line and 
the adjacent continuum, respectively. A five-point sub-pixel 
dither pattern (WFC3-UVIS-DITHER-BOX) was executed 
between individual exposures to improve sampling of the 
WFC3 point spread function (PSF). Since NaStl exhibits a 
small angular size, images were read out in a IKxlK sub¬ 
array mode to decrease readout overhead. To mitigate the 
effect of degrading charge transfer efficiency (CTE) in the 
WFC3 CCD chip, a 12 e~ post-flash was performed at the 
end of each exposure, which slightly elevated the background 
counts and mitigated against CTE losses. 

Some of the F658N and F645N images obtained on 
March 20 were severely degraded because the fine-guidance 
sensors failed to track the guide stars from the Guide Star 
Catalog, which turned out to be blends of multiple stars at 
HST resolution. The observations thus defaulted to GYRO 
pointing control and guiding, with no fine guidance. The 
resulting stellar PSFs, which normally have optimal full- 
width at half-maximum (FWHM) values of ~ 0707, were 

^ Note, however, that sensitivity of interferometry to spatial fre¬ 
quencies depends on the distribution of the baselines. 
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Figure 1. HST /WFC3 narrowband images of NaStl through F658N and F645N filters (logarithmic stretch), corresponding to the [N ll] 
emission line and adjacent continuum, respectively (upper panels). The lower panels show the F658N image with the central point source 
removed, and a 3-pixel smoothed version. 


highly elliptical and, in some cases, double-peaked, exhibit¬ 
ing FWHM values ranging from O'.'lO to 0”l9. Although 
several of the individual images from this visit were salvage¬ 
able and scientifically useful, the overall telescope drift dur¬ 
ing the exposures significantly limited our planned imaging 
sensitivity and resolution. NaStl was thus re-observed on 27 
April 2013 (UT 21:04), using an observing strategy identi¬ 
cal to the 20 March observations, although with new guide 
stars that were confirmed to be single. The repeated obser¬ 
vations achieved the intended sub-pixel dithering accuracy 
and yielded images with FWHM PSF cores of 0('09 using 
the Multidrizzle pipeline. Three of the lesser quality F658N 
[N ll] images and two of the F654N images from the initial 
observations were of acceptable quality and worth adding 
into our final drizzled image. 

Figure[T] shows the drizzled F658N and F654N images 
of NaStl. Two sets of diffraction spikes with different in¬ 
tensities can be seen, as a result of combining images from 
two separate visits with HST at different roll angles for the 
spacecraft. The nebula is clearly discernible in the F658N 
image, while no extended structure is immediately appar¬ 
ent in the F654N image. Subtraction of the bright cen¬ 
tral point source was challenging, and the most optimal 
strategy for getting the best possible subtraction was dif¬ 
ficult to converge upon. The F658N line and F654N con¬ 
tinuum filters produce PSFs that have substantially differ¬ 
ent speckle distributions and Airy-ring patterns, so simple 
subtraction of the continuum image from the line image 
did not produce satisfactory results, leaving large residu¬ 
als with complex mottled structure. Furthermore, both fil¬ 
ters exhibit unique ring-like ghost patterns that extend to 
« 0”9 from the PSF centroicQ Subtraction of artificial PSFs 

^ WFC3 Instrument Handbook (§6.5.3). 


generated by the Tiny Tim software also produced unsat¬ 
isfactory results, since the PSF models do not account for 
filter-dependent ghost artifacts. We thus found that the best 
strategy for subtracting the central point source was to use 
a bright neighboring star from the same images, 15" away 
from NaStl. We cut out a sub-image centered on this star, 
smoothed the pixels outside of the PSF halo to mitigate 
the effect of added noise to the subtracted image, and then 
scaled the ‘model’ incrementally to find the optimal subtrac¬ 
tion that minimized the difference residuals. In the end, the 
best scale factor of 0.8 removed all the higher-order resid¬ 
ual features in the subtracted image very well (the Airy 
rings, speckles, and diffraction spikes), but left a substantial 
amount of point-like [N ll] emission at the core of the neb¬ 
ula. Some of this point-like emission could be real, or could 
also be an artifact of imperfect subtraction. To check for 
the presence of compact continuum emission, we performed 
the same strategy for the F645N image and determined that 
there is no detection of an extended component from con¬ 
tinuum emission or scattered light at the image resolution, 
just noise residuals. 

Equatorial geometry 

The lower panels of Figure[T]show the star-subtracted F658N 
image and a 3-pixel smoothed version, which shows the inner 
structure of the nebula more clearly. Overall, NaStl’s nebula 
exhibits a complex multi-component structure having disk¬ 
like morphology with approximate major and minor axes of 
7" and 2", respectively, corresponding to an axis ratio «3.5. 
The emitting material thus appears to have equatorial dis¬ 
tribution and, if approximately circular, is tilted at an angle 
of « 12 ± 3'’ from edge-on (inclination angle of i = 78 ± 3°). 
The central point source is slightly offset by 60 mas from 
the geometric center of the nebula, as determined by the 
distribution of isophotal contours of the faintest outer re- 
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Figure 2. Left: star-subtracted 3-pixel smoothed F658N image of NaStl (logarithmic stretch), displayed in false color to enhance the 
nebular structure. Right: same as above with 9-pixel smooth and more extreme stretch to enhance the faint outer nebula. Flux contours 
of the brighter central regions from the 3-pixel smoothed image are overlaid for reference. The orientation of the images is the same as 
in Figure[T] 


gions. Figure[2] shows smoothed false-color versions of the 
star-subtracted [N ll] image. The inner portion of the nebula 
is composed of a relatively bright and smooth central ellip¬ 
soid with major and minor axes of 2'.'3 and 0”6, respectively, 
consistent with the same axis ratio as the greater structure. 
Based on the distribution of isophotal contours, the central 
point source also appears to be slightly offset from the center 
of this inner ellipse by 50 mas. Outside of the central ellipse 
there is a gap in length along the major axis, most 
notably on the south-west side. Beyond this gap, the outer 
nebula is composed of a clumpy ring or arm that can possi¬ 
bly be traced around the entirety of the structure. There is a 
particularly bright cluster or knot of emission clumps within 
this outer structure, at «l"4 south/south-west of the cen¬ 
tral point source. The brightest regions of the outer ring on 
either side of the center appear to have approximate point- 
reflection symmetry, which hints at a spiral structure for the 
outer ring/arm(s). The offset of the central point source from 
the apparent geometric center of the nebula could also be 
an indication of a one-armed spiral morphology. A highly 
contrast-stretched 9-pixel smoothed image (Figure[2j right 
panel) shows faint nebular structure with a full diameter of 
Y!7. 


2.2 Optical spectroscopy 

Optical spectroscopy of NaStl was performed on 2012 March 
31 using the 2.3 m Bok Telescope at Kitt Peak, Arizona. A 
low-resolution grating having 4001mm“^ was used in con¬ 
junction with a T.'5 slit to achieve a FWHM spectral resolu¬ 
tion of 4.6 A and a wavelength range of AA4150-7200 A. Four 
exposures of 120 s each were obtained at an airmass of 1.5- 


1.6 and median combined into a final spectrum. Flat fielding 
and wavelength calibration were performed by taking spec¬ 
tral images of the internal continuum and arc lamps at the 
beginning of the night. Flux calibration was performed us¬ 
ing spectra of the white dwarf standard star BD+42 2811. 
Spectroscopy was not performed with the slit aligned to the 
parallactic angle, so some loss of blue light is expected. All 
image and spectral reductions were conducted using stan¬ 
dard IRAF routines. 

Spectroscopy of NaStl was also performed on 2012 Oct 
17 with the 6.5 m telescope at MMT Observatory using the 
Bluechannel spectrograph. The 1200 1 mm“^ grating was uti¬ 
lized with a 1” slit to achieve a dispersion of «0.5 A per pixel 
(«20kms“^ pixel”^). For this observation, the slit was ori¬ 
ented at 31° from north, approximately parallel to the major 
axis of the nebula. This strategy was used in order to search 
for kinematic features that might be discernible along the 
resolved axis of the nebula. Three exposures of 300 s each 
were obtained at an airmass of 1.4 and median combined 
into a final spectrum. Flat helding and wavelength calibra¬ 
tion were performed by taking spectral images of the internal 
continuum and arc lamps at the beginning and end of the 
night. Flux calibration was performed using spectra of the 
white dwarf standard star Feige34. Since the slit was not 
aligned with the parallactic angle, some loss of blue light is 
expected. 

Figure|3] shows the low- and moderate-resolution optical 
spectra of NaStl from Bok and MMT. The spectral mor¬ 
phology is mostly nebular in nature, exhibiting a complex 
mix of permitted and forbidden emission lines spanning a 
broad range of ionization. The comparatively higher reso¬ 
lution spectrum from MMT shows the resolved NalD dou- 
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Figure 3. Optical spectra of NaStl from Bok (upper spectrum) and MMT (lower spectrum). The most prominent nebular features have 
been marked (see CS99 for a comprehensive list of all emission lines). 
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Figure 4. Radial velocity of the [N ii] A6584 profile (left panel, solid curve) as a function of position along the slit (middle panel), which 
was approximately aligned with the major axis of the nebula, and the same for an arc line of Ne 1 6599 A (right panel). Real kinematic 
motion for the outer nebula of NaStl is apparent. Note that velocity measurements for offsets < 2" are dominated by emission from the 
central point source, and are not reliable. 


blet from interstellar (and, perhaps, circumstellar) absorp¬ 
tion along the line of sight. The MMT spectrum also shows 
extended Lorentzian wings for the Ha, close examination of 
which shows a full possible extent of at least ±2000 km s~^ 
for the wing edges; this Lorentzian profile indicates that Ha 
is not entirely nebular in origin, and suggests that some of 
the flux is from a physically distinct region of high optical 


depth, where electron scattering could be occurring. Nar¬ 
rower wings might also be associated with the He A6680 A 
profile, but not at the level observed for Ha. Overall, the op¬ 
tical spectra presented here are consistent with the previous 
observations of CS99, and we find no substantial spectral dif¬ 
ferences between the epochs presented here or with respect 
to the measurements in the literature. 
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Figure 5. Top-, spatial profile of selected emission lines along the 
spectrograph slit, which was tilted near the long axis of the nebula 
at 35°. The lines are scaled in intensity to roughly match the peak 
of the [N ll] A6584 feature. The nebula is clearly resolved in the 
[N ll] transitions. The bumps on either side of the continuum 
profile are an artifact of the scaling up of the continuum, and are 
not real resolved continuum structure. 


Spatially resolved kinematics 

Figure|4] (center panel) shows a two-dimensional spectral im¬ 
age of the [N ll] A6584 line from the MMT observation, for 
which the slit was aligned along the major axis of the neb¬ 
ula. A similar image of a neighboring arc line (right panel) 
from the calibration exposures is also shown for comparison. 
The left panel of Figure|4] shows the measured shift in ra¬ 
dial velocity of [N ll] along the spatial axis, which exhibits 
an ‘S’ shaped curve having a peak velocity of ±10 kms“^ 
at from the center of the continuum dispersion axis 
(we note that the velocity measurements for offsets < 2" 
are dominated by emission from the PSF of the central 
point-source, and are thus not reliable tracers of the neb¬ 
ular velocity at those positions). We do not see any such 
velocity structure associated with Ha or Hel; this is to be 
expected, since these features are spatially unresolved. For 
example, Figure[5] shows the spatial profile of [N ll] along the 
slit compared to Ha, Hel, and the continuum adjacent to 
[N ll] A6584. It is apparent that the [N ll] emission from the 
nebula is resolved, while the other prominent lines and the 
continuum are not. 

It thus appears likely that the [N ll] emission exhibits 
real kinematic structure that is intrinsic to the object, likely 
to be an imprint of outflowing gas. If the nebula is an ex¬ 
panding disk or ring then the relative motions of gas along 
the slit, which are aligned to the major axis of the neb¬ 
ula, will be mostly transverse, not radial. In this case, the 
true expansion velocity could be significantly higher. Spec¬ 
troscopy with high spatial and spectral resolution, sampled 
along the minor axis of the nebula is necessary to gauge the 
true velocity structure. 


3 INFRARED IMAGING AND 

SPEGTROSCOPY 

3.1 AO imaging with MagAO±Glio2 

Near-infrared AO images were obtained on 2013 April 7 with 
the Magellan 6.5 m Clay telescope at Las Campanas Ob¬ 
servatory during commissioning of the MagAO and Clio2 
instrument (Morzinski et al. 2014). The wide camera was 
used, yielding a plate scale of 27.5 mas and a field of view 
of 28” X 14”. NaStl itself served as the natural guide star 
for the AO system. Thirty images of 40 coadds each were 
obtained through a Ka filter and 10 images of 40 coadds 
each were obtained through L' and M' filters. The L' and 
M' observations also included a neutral-density (ND) fil¬ 
ter to prevent saturation of the detector. We note that 
the Ks filter setup used throughout the observing run did 
not have a long-wavelength blocker, which resulted in a 
small amount of leakage of light from 3-5 /rm wavelengths 
(hereafter, reference to the Ka band includes the potential 
for some long-wavelength leakage). The telescope was nod¬ 
ded by 4” between exposures for background subtraction. 
No fiat-fielding was performed. The background-subtracted 
data were then median combined into final image frames for 
each filter. Three field stars in the Clio2 image have entries 
in the United Kingdom Infrared Deep Sky Survey (UKIDSS; 
Lawrence et al. 2007) Galactic Plane Survey (GPS) catalog 
(Data Release 7), and were used to calibrate the astrome¬ 
try. Two of those stars were also present in the HST /WFC3 
frame, which enabled us to align the astrometry of the op¬ 
tical and infrared images to a relative accuracy of 0V064 in 
physical arcsec. 

Figurein] (upper three panels) shows the MagAO/Glio2 
images of NaStl in the three infrared bands. At M' NaStl 
appears to be unresolved, as the diffraction-limited structure 
of the circular PSF is clearly discernible in the form of mul¬ 
tiple concentric Airy rings. The FWHM of the M' PSF cen¬ 
troid is 5.96 pixels, which corresponds to an angular size of 
164 mas. At L' the FWHM is 142 mas and the first Airy ring 
is marginally discernible. The PSF halo is asymmetrically 
bright, which could indicate that the source is marginally 
resolved. The As-band image, by comparison, shows that 
the source is fully resolved into an ellipsoidal structure that 
is aligned with the halo asymmetry observed at L'. The 
FWHM of the Ka centroid along the longest axis is 171 mas, 
substantially broader than the PSF centroids of the L' and 
M' images. It thus seems unlikely that the ellipsoidal struc¬ 
ture in the Ka band has resulted from the small amount of 
potential leakage from 3-5 /rm light. Moreover, observations 
of an unresolved star (2MASS J12073100—3932281) having 
the same Strehl ratio produced a substantially more com¬ 
pact stellar Ka PSF having FWHM=123mas. In addition, 
observations of HD 73634, which is comparably bright to 
NaStl in the 3-5 /rm range, produced a stellar PSF having 
FWHM=126mas. Any potential degradation of the stellar 
PSF from long-wavelength leakage, if present, should have 
been noticeable in the images of this star as well. No such 
effect is observed, and we thus conclude that the extended 
As-band emission around NaStl is real. At L' and M', how¬ 
ever, NaStl has the same FWHM size as the comparison 
stars and thus appears to be unresolved at these longer 
wavelengths, although it is possible that there is marginally 
resolved structure in the halo of the L' source. 
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Figure 6. Upper panels: MagAO Ks, L', and M' images of NaStl 
in false color. The images are ^2'.'7 square in angular size and ori¬ 
ented in detector coordinates. Lower panels: the ifs-band image 
of NaStl with logarithmic intensity contours (left) and the HST 
[N ll] narrow-band image from Figure 2 with the Ks contours 
overlaid (right). An ellipsoid of extended Kg emission is clearly 
resolved, and is approximately aligned with the major axis of the 
optical nebula. The lower images are oriented with north up and 
east toward the left. 


FigureE] (lower panels) also compares the extended Ks 
image of NaStl to the [N ll] image from HST. Interestingly, 
the position angle of the Ks ellipsoid on the sky is roughly 
parallel to that of the major axis of the optical nebula, 
which provides further reassurance that the extended in¬ 
frared structure is real. 

The FWHM of the long axis in the Ks band (171 mas) 
corresponds to a physical diameter of 171-513 AU, given the 
range of possible distance (1-3 kpc), and the outmost regions 
of the Ks nebula could extend to 700-2100 AU. The near- 
infrared angular size is thus substantially larger than the 
8-13 pan size measured from the interferometric imaging of 
Rajagopal et al. (2007). Those authors found that the an¬ 
gular size of NaStl increases from 15 to 20 mas between 8 
and 13 pm, respectively. This either implies that their inter¬ 
ferometric measurements were not sensitive to lower spatial 
frequencies from relatively large angular sizes of >140 mas, 
or that the apparent size is of NaStl is intrinsically larger 
at 2-3 pm than at 8-13 pm, and perhaps not the result of 
the same emission source. 


3.2 SINFONI spectrum 

We obtained infrared spectra of NaStl from the European 
Southern Observatory (ESO) archive. The object was ob¬ 
served with the ESO VLT (Yepun) telescope on 2005 Aug. 
30 using the integral field spectrograph SPIFFI within the 
SINFONI module (Eisenhauer et al. 2003). The instrument 
delivers a simultaneous three-dimensional data cube with 
two spatial dimensions and one spectral dimension. The K- 
band (1.95-2.45 pm) grating was used with a resolving power 
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Figure 7. A'-band spectrum of NaStl from SINFONI. 


R ~ 4000. The pixel scale was d'.l with a field of view 
of 3 V 2 square. Observations were performed under seeing- 
limited conditions, which were relatively poor at 1.4-1.6 arc- 
sec. The airmass during the observations was 1.1. NaStl was 
thus spatially unresolved by SPIFFI. The data were reduced 
using the publicly available ESO SINEONI pipeline. The 
A'-band spectrum was extracted by averaging all “spaxels” 
within an aperture that encompassed the «l'.'5 seeing disk 
of the unresolved source. Telluric correction was performed 
using a SINFONI spectrum of the B5 star HD 183734 from 
nearly the same airmass. 

Figure!?] shows the A'-band spectrum of NaStl. The 
morphology and relative line strengths appear consistent 
with the spectrum presented in CS99. Overall, the spectrum 
is dominated by an extremely strong emission line of He l at 
A2.06 pm, and the line profile exhibits Lorentzian wings ex¬ 
tending to ~ ±2500 km s”*^, similar to the Ha line in the 
MMT spectrum (see FigureEJ. This morphology suggests 
that some of the emission is emerging from a region of high 
optical-depth, where electron scattering is occurring. The 
total integrated flux across the line amounts to «12% of the 
integrated counts between 2.0 and 2.3 pm (before telluric 
correction). The spectrum also exhibits strong He l emission 
at 2.113pm, a blend of Bry and Hel near A2.166pm. Hell 
at A2.189pm, and Fell at A2.089pm. 


3.2.1 The nature of the extended Ks-band emission 

Based on the infrared SED presented in CS99, NaStl is 
clearly a strong source of thermal infrared excess emission 
from dust, which is composed of at least two dust compo¬ 
nents having respective temperatures of 700 K and 2000 K. 
It is thus natural to consider this dust as a potential source 
of the extended flux resolved by MagAO in the Ks band. 
However, if the extended flux is thermal continuum emission 
from dust, then the source should have also been resolved 
at M' (where NaStl is brightest) for any reasonable range 
of dust temperatures. Indeed, at >171 AU distances from 
the central heat source, hot dust formed within an outflow 
should have a cooler M' counterpart with a larger angular 
size. Instead, the PSF centroid at M' is substantially smaller 
than the angular size of the resolved source in the Ks band, 
and no such extended counterpart is detected. This raises 
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doubts that the extended Ks component results from hot 
dust emission. 

There is an alternative possibility that the extended 
ffs-band emission is the result of photoionization. Given the 
extremely bright flux of the Hel emission lines observed in 
the if-band spectrum, and the approximate sky alignment 
between the Ka and [N ll] extended structures, it is plausi¬ 
ble that the MagAO observations resolved the nebula in the 
light of He I, the flux of which is dominated by the A2.06 /rm 
transition. This also suggests that the marginally resolved 
structure in the L' PSF halo could be extended Hel emis¬ 
sion at A3.703/rm. The compactness of Hel relative to [N ll] 
would explain why NaStl was not spatially resolved from the 
ground in the light of He l by our spectroscopy (see Figure[S|, 
nor by the narrowband imaging of CS99. However, the broad 
Lorentzian wings in the Hel A2.06/im line prohle indicate 
that electron scattering is occurring, which is not likely to be 
efficient in the extended nebula at >171 AU distances from 
the central star. Rather, this type of line profile is sugges¬ 
tive of emission from a dense stellar wind, which would have 
to be produced near the central star and, hence, in the un¬ 
resolved core of the nebula. Thus, if the broad-winged Hei 
emission is responsible for the resolved flux at Ks, then it 
could be emission from the central stellar wind further scat¬ 
tered into the outer resolved regions by dust that is too cool 
to emit its own detectable continuum in the near-infrared. 
This could explain the larger angular size at Ks versus M', 
since dust scattering is more efficient in the former band, as 
has been observed in rj Car’s Homunculus nebula (Smith et 
al. 1998). In this case, however, it is puzzling that no optical 
counterpart of dust-scattered fight was also detected in our 
HST F654N optical continuum image of NaStl. 

Unfortunately, we cannot attribute the Hel A2.06pm 
line profile to specific regions within the nebula resolved by 
MagAO, since the source was spatially unresolved by the 
SINFONI observations. The broad-winged profile could sim¬ 
ply be from the inner unresolved portion of the nebula, and 
a narrower line profile could comprise the extended nebula 
resolved at Ks. Future AO-assisted spectroscopy and near- 
infrared spectropolarimetry could help discriminate between 
direct line emission from the resolved nebula and dust scat¬ 
tering of stellar wind emission emerging from the unresolved 
interior. 

We conclude that the extended flux at Ks is most likely 
to be fight from Hel A2.06pm, where possibly some, and 
perhaps most, of the flux could have been scattered out to 
>171 AU by dust too cool to emit strongly at Ks, L', and 
M' . Therefore, the strong point-source emission at these 
wavelengths must be from a compact component of hot dust 
near the core, perhaps physically associated with the mid- 
infrared source resolved by interferometry (Rajagopal et al. 
2007). 


4 CHANDRA X-RAY OBSERVATIONS 

We observed NaStl on 2014 September 24 19:01 UT using 
the Chandra X-ray Observatory (Obs. ID 15700) and the 
Advanced CCD Imaging Spectrometer (ACIS) instrument 
(Weisskopf et al. 2002). A total exposure time of 12 kilosec- 
onds resulted in the detection of a faint point source at a = 
18:52:17.53 5 = -f00:59:44.6 (J2000.0), with an positional 


error of l'.'25 (90% confidence), consistent with the respec¬ 
tive optical and infrared positions of the bright central point 
source from HST and Magellan. The proper designation for 
the X-ray source is thus CXOU J185217.53-1-0059446, in ac¬ 
cordance with the standard naming convention for Chandra 
sources. 

We extracted the source photometry and spectra using 
the CIAO routines SRCFLUX and SPECEXTRACT, respectively. 
After accounting for the background rate, the total num¬ 
ber of counts attributable to the object is 42 ± 7 in the 
broad 0.5-7.0 keV energy band. The associated photon flux 
and energy fluxes in this band are 1.2)';® x 10“®cm“®s“^ 
and 2.93;7 x 10“^“^ ergcm“®s“^, respectively (values accom¬ 
panied with 90% confidence envelopes defined by the sub¬ 
script and superscript), and the average energy of the de¬ 
tected emission is (E) — 1.7 ± 0.3 keV. The fluxes within 
the individual soft (0.5-1.2 keV), medium (1.2-2.0 keV), 
and hard (2.0-7.0 keV) bands are listed in Tablet, along 
with the hardness ratios, defined as HR^sd-soft = {Fmed — 

Faoft)/(Fsoft -t-Fmed -fFiard) and similarly for HRhard-med- 

The hard band exhibits the largest energy flux, while the 
majority of counts are in the medium band between 1 and 
2 keV. 

The ACIS-I energy spectrum of NaStl is shown in 
Figure|8] along with an ACIS spectrum of the WN4 star 
WR18 for comparison (vertically scaled by 0.5), which we 
obtained from the Chandra Data Archive (Obs. ID 8910, PI: 
S. Skinner). Overall, the spectra of NaStl and WR18 appear 
similar, both having the majority of counts between 1 and 
2 keV, heavy absorption at energies < 1 keV, and compara¬ 
ble flux at harder energies relative to their 1-2 keV fluxes. 
WR18 exhibits strong Sixill triplet emission at 1.84-1.86 
keV; this is a common feature in the X-ray spectra of WN 
stars that indicates emission from a thermal plasma hav¬ 
ing a maximum fine power temperature of logTmax = 7.0 K. 
NaStl also appears to exhibit an excess of flux at this energy; 
we suspect this may result from the same Sixv transition, 
although the detection is of relatively low statistical signifi¬ 
cance. WR 18 also exhibits a Sixv feature at 2.46keV, which 
indicates the presence of an additional hot plasma compo¬ 
nent having a higher temperature (logT^aa; = 7.2 K); this 
line is not detected in NaStl, although our relatively weak 
observational sensitivity at these higher energies does not 
allow us to place meaningful constraints on its presence. 

Skinner et al. (2010) modelled the full spectrum of 
WR 18 with a two-temperature thermal plasma having 
kTcooi = 0.65 ± 0.22 keV, kThot = 3.5 ± 1.9 keV, and an 
X-ray luminosity of Lx « 3 x 10®® ergcm“® s“^. Based on 
the observed similarities, it is plausible that the X-ray emis¬ 
sion from NaStl is also emitted by a thermal plasma having 
physical parameters similar to WR 18. Unfortunately, we are 
unable to attribute a unique plasma model to NaStl with 
any confidence because of the low number of counts (e.g., 
a power-law model can be tuned to fit the spectrum with 
the same level of statistical significance as a thermal plasma 
model, such as that used in the case of WR18). We can 
nonetheless place some constraints on the parameters of the 
emitting plasma using the broad-band X-ray photometry 
and assuming a plausible range of emission models. First, to 
estimate NaStl’s absorption-corrected X-ray flux, Fu„aba, 
we adopted the extinction value of Ay ~ 6.5 mag from 
CS99 to estimate the total hydrogen absorption Nh column 
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Figure 8. Chandra ACIS-I imaging spectra of NaStl and the 
WN4 star WR18 (red, Skinner et al. 2010) scaled by a factor 
of 0.5. The excess of flux for at 1.8—l.QkeV in the spectrum of 
WR18 is due to the Sixill triplet; NaStl exhibits a marginally 
significant excess at the same energy. 


along the line of sight via the Av/Nh relation from Pre- 
dehl & Schmitt (1995), and obtained Nh ~ 1.2 x 10^^ cm^. 
We applied this value to several optically thin thermal 
plasma models having single temperatures in the range of 
kT = 0.6-2 keV (a reasonable range for WRs), and obtained 
Funabs = (7-38) X10“^"* ergcm“^ s”*^. We also applied a sim¬ 
ple power-law model having a photon index in the range 
of r=l-2, and obtained a smaller range of Funabs = (3- 
6 ) X10”^“^ erg cm~^ s“^. Given the likely presence of a WR 
in NaStl and the observed similarities with the X-ray prop¬ 
erties of WR 18, we favor a thermal plasma model for this 
source. If the flux we detect is mostly from a 0.65 keV 
plasma component, like that of WR 18, then the unab¬ 
sorbed flux is ~ 7 X lO”*^^ ergcm“^ s“^. For a distance in 
the range of 1-3 kpc, the corresponding X-ray luminosity 
would be Lx ~ 10®^-10®® ergcm“^ s“^. Given the poten¬ 
tial for additional self absorption to X-rays, this should be 
regarded as a plausible lower limit to the X-ray luminosity. 
Adopting the stellar bolometric luminosity range from GS99 
(4 < logL/L© < 6.5), the X-ray luminosity estimate implies 
Lx/Lboi 10 “®- 10 “®. 

To consider the full possible extreme for self-absorption 
we could compare NaStl to the subset of WR binaries 
whose stellar spectra have also been completely veiled by 
an optically thick dust screen (e.g. the DWGL “cocoon” 
stars that populate the Quintuplet Cluster in the Galac¬ 
tic Center; Figer et al. 1999; Mauerhan et al. 2010). Such 
objects are probably near the high end of self-absorption for 
massive stars. For example, CXOGC J174645.2—281547, lo¬ 
cated in the Galactic center field, is consistent with an op¬ 
tically thin thermal plasma having a high temperature of 
kT « 3.8 keV heavily absorbed by an intrinsic column of 
Ah = 2 X 10^^ end (Hyodo et al. 2008). If we apply such 


Table 1. Basic X-ray parameters for NaStl 


Name 

eXO J18521753-I-0059446 

Fbroad(0.5-7.0 keV) 

2.9^'^ X lO-i'^ 

F,oft(0.5-1.2 keV) 

4.8| f X 10-1® 

Fn,ed(1.2-2.0 keV) 

l.llf X lO-i-i 

Fhard(2.0-7.0 keV) 

CO 

oto 

X 

o 

1 

^ ^med — soft 

0.15 ±0.14 

^ ^hard — med 

0.18 ±0.17 

(E) 

1.8±0.3keV 

Ah (ISM) 

1.2 X 10^2 cm2 

-^X 

~ 10®2-1033ergs-l 

Lx/Lbol 

~ 10“®-10“® 


Notes. Energy fluxes are in nnits of ergcni~^s~^ and the snperscripts 
and subscripts represent the 90% confidence envelopes. The hardness 
ratios so/t ) were computed using photon fluxes to calculate 

(^’med - Aoft)/(Aoft + turned + t^hard) and similarly tor ■ 

The interstellar hydrogen column Nj^(ISM) was computed by adopting 
Ay = 6.5 mag (CS99) and using the relation of Predehl Sc Schmitt (1995). 
X-ray luminosity Lx was calculated using distance range of 1—3 kpc and 
assuming a range of intrinsic X-ray spectra properties (see text §3). 

a value to NaStl and assume the X-rays are also produced 
by optically thin thermal plasma having a temperature in 
the broad range of 0.6-3.8 keV, then we obtain unabsorbed 
fluxes in the range of (1-100) xl0“^® ergem”^ s“^. The cor¬ 
responding range of possible luminosities would therefore be 
Lx 10®^-10®^ ergs”^. However, it appears that the intrin¬ 
sic plasma temperature of the dominant emission compo¬ 
nent in NaStl is unlikely to be higher than 2 keV, given 
the observed drop in X-ray flux at energies above this value 
and the overall similarity to the cooler WR 18. Furthermore, 
an absorption column as high as Ah = 2 x 10 ^® cm^ ap¬ 
pears to be inconsistent with the data as well, since all 
of the emission models we have considered (optically thin 
thermal plasma, thermal Bremsstrahlung, and power law) 
would have their emission strongly suppressed at energies in 
the range 1-2 keV if the absorption were this strong. Most 
of the X-ray flux from NaStl is actually in the range of 
1-2 keV, so it seems unlikely that the absorption column 
is greater than 10^^ cm^. It thus appears that a total col¬ 
umn of Ah = 2 X 10^^ cm^ inferred from the adopted total 
optical-infrared extinction might be a close approximation, 
and could already account for a substantial fraction of the 
self absorption. In consideration of all of the facts, we con¬ 
clude that the X-ray luminosity of NaStl is most likely to 
be Lx ~ lO^^-lO^® ergs“^; but given the potential for the 
presence of a fainter hard component that the observations 
were not fully sensitive to, and the potential for some addi¬ 
tional self absorption, the actual value is probably closer to 
lO^^ergs”^. In turn, this implies Lx/Lboi ~ 10“^-10“®. 

The X-ray properties of NaStl are broadly consistent 
with the X-ray characteristics of both purportedly single 
WRs and those in binaries (Skinner et al. 2010; Gagne et al. 
2012). The hot thermal plasma in these systems is generated 
by shocks that result from supersonic winds, in many cases 
by two opposing winds in collision. WR 18, which we com¬ 
pared to NaStl in §4, is not a confirmed binary star, so it 
not entirely clear if the X-ray properties in that system or in 
NaStl are the result of colliding winds. Indeed, many other 
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colliding-wind binaries (CWBs) have more prominent hard 
X-ray components {kT > 2keV, Gagne et al. 2012), and 
in some cases exhibit Fe Ka emission near 6.7 keV energy, 
for example, in rj Car (Corcoran et al. 2000) and WR 147 
(Skinner et al. 2007). As noted earlier, the lack of a simi¬ 
larly hard component in NaStl could simply be due to the 
current observational sensitivity limits. Intrinsic dimness at 
hard X-ray energies, on the other hand, could indicate that 
the winds of the component stars are relatively slow com¬ 
pared to the stars in more hard-X-ray bright CWBs (Luo 
et al. 1990; Usov 1992). Hard X-ray emission can also vary 
with orbital phase, so NaStl could have been observed in a 
relatively soft state. Deeper X-ray observations in the future 
will allow for more meaningful constraint on the presence of 
hard X-rays and variability in NaStl. 


5 DISCUSSION 

CS99 suggested that NaStl is a WR cloaked in an opaque 
dusty nebula, similar to the Homunculus around r) Car. 
Indeed, their nebular analysis demonstrated that NaStl is 
the only other object in the Galaxy known to exhibit such 
an extreme level of nitrogen enrichment (indicative of ad¬ 
vanced CNO processing from an evolved massive star), while 
their derived parameters for the ionization source imply 
T >> 30 kK and log (L/Lq)=5.1-6.1, consistent with a WR 
star. Our discovery of X-ray emission from NaStl, shown in 
§4, provides additional evidence for a WR, where hot plasma 
is produced by shocks in supersonic stellar winds, perhaps 
from colliding winds in a binary. However, although NaStl 
shares some interesting similarities with ry Car’s Homuncu¬ 
lus, the new data highlight some obvious differences, de¬ 
scribed as follows. 

The Homunculus was created during an extremely ener¬ 
getic explosion that ejected at least « 12 Mq of material at 
velocities of « 600kms“^ (a ~ 10®° erg event; Smith et al. 
2003). The outflow exhibits a striking bipolar morphology 
and is a strong source of [N ll] emission, but also Ha and 
strong optical-infrared continuum from scattered light and 
thermal radiation from dust. The extended nebula around 
NaStl, on the other hand, is clearly an equatorial structure, 
with no obvious sign of a bipolar component. The optical 
emission is pure [Nil], with no extended Ha emission or 
extended continuum evident in our HST images (Figure 1) 
or ground-based spectroscopy (Figured). Although we have 
yet to obtain spatially-resolved kinematic spectroscopy of 
NaStl along the minor axis of the nebula, where we ex¬ 
pect the outflow’s apparent radial velocities to be highest, 
there is currently no indication that the outflow velocities 
are anywhere near those of rj Car’s explosive Homunculus. 
Unfortunately, we are unable to obtain a reliable mass esti¬ 
mate for NaStl’s nebula based on [N ll] alone, so information 
on the kinetic energy of the outflow is lacking. Nonetheless, 
the physical processes that created NaStl’s nebula appear 
to have been far less energetic than the giant eruption that 
created rj Car’s Homunculus, based on the available data. 
Below we speculate on the origin of NaStl in the context of 
other known massive binaries that have equatorial nebulae 
and exhibit X-ray emission and thermal dust emission. 


5.1 Dust formation from colliding winds? 

NaStl’s equatorial nebula and X-ray-emitting plasma, and 
the compact thermal infrared component detected by CS99 
and resolved by Rajagopal et al. (2007), are reminiscent 
of a subset of WR binaries that produce dust (Gehrz & 
Hackwell 1974; Williams et al. 1990; Oskinova & Hamann 
2008), episodically in some cases (Williams et al. 1978, 1990; 
Veen et al. 1998). The ring-like and possible spiral structure 
in NaStl brings to mind the specific cases of WR98a and 
WR 104 (Monnier et al. 1999, 2002; Tuthill et al. 1999, 2008), 
which are CWBs containing WC stars surrounded by strik¬ 
ing Archimedean spirals of hot dust emission viewed nearly 
face on (the so-called ‘pinwheel’ nebulae). The dust in these 
systems is thought to be produced when a carbon-rich wind 
from a WC star mixes with the H-rich wind of a less-evolved 
OB companion, creating chemistry suitable for the forma¬ 
tion of graphite or silicates. After strong shock compression, 
dust condenses within a conic cooling flow that develops 
downstream from the wind interaction zone (Williams 2008). 
It is interesting in this regard that the thermal infrared emis¬ 
sion comprising the spiral nebulae of WR 104 and WR 98a 
have physical diameters of a few x 100 to 1000 AU, because 
this is comparable to the infrared size we have measured for 
NaStl. It thus appears plausible that the resolved structure 
of infrared emission revealed by our AO imaging could also 
be a spiral structure viewed 12 ° from edge-on, whether the 
emission is from thermal dust continuum or scattered line 
emission. Meanwhile, the very recent formation of the hot 
dust, as suggested by Rajagopal et al. (2007), suggests the 
possibility of new evolutionary changes that only recently 
created suitable conditions for dust formation. One possi¬ 
bility is that the central WR recently transitioned from the 
WN to the WC phase after the ejection of N-rich material. 
Indeed, there is a growing body of evidence from observa¬ 
tions of SNe Ibn, which indicates that WN stars can undergo 
discrete super-Eddington eruptions before exploding (Foley 
et al. 2007; Pastorello et al. 2008). If this is the case for 
NaStl, the eruption could have been followed by the devel¬ 
opment of a C-rich wind from the newly formed WC star. 

C-rich chemistry is not a requirement for dust to form 
in a CWB. The circumstellar material around ij Car, for ex¬ 
ample, which is N-rich yet C-poor, exhibits firm evidence 
for freshly synthesized hot dust {T = 1400-1700 K). Sim¬ 
ilar to WC binaries, infrared monitoring has shown that rj 
Car forms dust episodically, apparently maximized by strong 
shock compression of the stellar wind during periastron en¬ 
counters (Smith 2010). In this case, the post-shock cool¬ 
ing apparently results in the formation of corundum dust 
(AI 2 O 3 ), rather than graphite or silicates. However, corun¬ 
dum usually produces a spectral emission feature near 12 jim 
wavelength, which is not seen in the relatively featureless 
spectrum of NaStl presented by Smith & Houck (2001). Sil¬ 
icates are also expected to produce spectral features in the 
8-14 jim range, unless the grains are extremely large. Amor¬ 
phous carbon, on the other hand, can produce a feature¬ 
less mid-infrared continuum (McDonald et al. 2012). Amor¬ 
phous carbon has also been favored in the cases of other dust 
producing WRs (Williams et al. 1987), and so it might be 
the most plausible dust chemistry for NaStl. However, the 
limited 8-14 ^m wavelength coverage of the existing mid- 
infrared spectrum could make very broad spectral features 
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from various dust species difficult to identify, and it might 
be premature to characterize the mid-infrared SED as “fea¬ 
tureless” based on those results alone. Broader wavelength 
coverage of the mid-infrared spectrum will thus be necessary 
to elucidate the dust chemistry in NaStl. 

Another open question is whether NaStl forms dust 
episodically, similar to other WR CWBs (e.g, see Williams 
et al. 1992, 1994) and 77 Car (Smith 2010). Such behavior 
would indicate an eccentric orbit for the binary. We com¬ 
pared the 12 fim flux of NaStl from the Wide-field Infrared 
Survey Explorer { WISE) AZZuiise catalog (Wright et al. 2010) 
with the flux measured at the same wavelength by Smith 
& Houck (2001) over a decade earlier; no change in mid- 
infrared brightness was detected. However, this does not 
rule out episodic dust formation, which is usually evidenced 
by variability in the near-infrared, since the freshly synthe¬ 
sized dust is relatively hot. As such, mid-infrared variations 
are not necessarily expected, since emission at these longer 
wavelengths will come from cooler dust that has flowed away 
from the system and piled up, to become more steadily illu¬ 
minated by the luminosity of the central stars. Interestingly, 
van der Hucht et al. (1997) noted near-infrared variability of 
NaStl at the level of 0.6 mag over a month-long timescale, so 
the possibility of episodic dust formation is worth investigat¬ 
ing with continued near-infrared monitoring in the future. 
Episodic dust formation might also be correlated with varia¬ 
tions in X-ray luminosity and hardness from colliding-wind 
emission (e.g., Williams et al. 2013), which is also worth 
exploring. Such data could reveal an orbital period for the 
system, which would provide much needed insights into the 
physical parameters of the component stars. 

5.2 Is NaStl the product of mass transfer? 

NaStl’s nebula of slowly moving CNO-processed mate¬ 
rial could also have become an equatorial structure as 
the result of binary mass transfer. If the evolved primary 
star recently underwent a large radial expansion, material 
could have been transferred to a companion via RLOF. In 
the case of non-conservative mass transfer, high angular- 
momentum material can bypass accretion and escape the 
system through the outer Lagrange points (Vanbeveren et 
al. 1998). This might explain the relatively low velocity of 
the outer nebular material (compared with WR winds) and 
its confinement to the equatorial plane. The Roche lobe need 
not be filled by the star itself. In relatively wide binaries 
mass transfer can also occur if the donor overhlls it Roche 
lobe with a slow dense wind (WRLOF; Abate et al. 2013), 
and the same is probably true for dense material from a dis¬ 
crete mass eruption, if the material moves slowly enough. 
Hydrodynamic simulations of relatively low-mass RLOF bi¬ 
naries containing asymptotic giant branch (AGB) stars have 
shown that a circumbinary envelope with spiral structures 
and diffuse clumpy rings can form in the equatorial plane 
and enshroud the system (Sytov et al. 2009). We speculate 
that NaStl could be a scaled-up analog of this process for a 
more massive binary. 

Mass transfer might also explain some of the kinematic 
features previously observed in NaStl’s spectra, in particu¬ 
lar, the double-peaked emission profiles detected by CS99. 
Those authors showed that the emission lines of various 
atomic species spanning a wide range of ionizations states 
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Figure 9. High-resolution normalized spectrum of the [N ll] pro¬ 
file of the central point source of NaStl, as a function of radial 
velocity (black solid line; reproduced from CS99), compared to 
that of the sgB[e]/LBV binaries Wdl-9 (red dashed line; repro¬ 
duced from Clark et al. 2013 ) and MWC314 (blue dash-dotted 
line; from Lobel et al. 2013), at arbitrary orbital phases, and RY 
Scuti (green dash-triple-dotted line; from Smith et al. 2002). The 
velocities were shifted to center the entire profiles. In the cases 
of Wdl-9 and MWC 314, the double peaked profiles have been 
interpreted as the result of rotational motion in a circumbinary 
disk, but could also be from an expanding ring or torus, as is the 
interpretation in the case of RY Scuti 

exhibit an interesting correlation between profile width and 
ionization potential, with the highest excitation lines ex¬ 
hibiting the largest velocities (FWHM up to ~160kms“^). 
This implies that there are multiple line-forming regions in 
NaStl’s circumstellar envelope. These characteristics pro¬ 
vide some important clues on the nature of NaStl, and fa¬ 
cilitate a comparison with other massive binaries in which 
mass transfer is thought to have recently occurred. 

Figure[9] shows the double-peaked [N ll] emission-line 
profile of NaStl (reproduced from CS99) and for several 
other binary systems containing evolved massive stars (note 
that the spectroscopic observations which produced this pro¬ 
file in NaStl were centered on the compact central source 
and did not sample the larger nebula). The separation be¬ 
tween each peak in the double profile is «30kms“^, compa¬ 
rable to other systems, and is particularly similar to the ob¬ 
ject Wdl-9 from the starburst cluster Westerlund 1 (Clark 
et al. 2013). Wdl-9 is composed of a supergiant B[e] star 
in a massive interacting binary that is undergoing, or has 
recently exited from, a phase of rapid RLOF. The star ex¬ 
hibits a rich emission-line spectrum with no photospheric 
features, similar to NaStl, albeit with softer ionizatior0. In 
addition to the [N ll] profile and overall spectral morphol- 

® Based on the presence of Hell emission, van der Hucht et al. 
(1997) suggested an 0[e] classification for NaStl 
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ogy, Wdl-9 is similar to NaStl in several other respects: it 
is an X-ray source having Lx ~ 10®®ergs“^ and plasma 
temperature consistent with colliding-wind shock emission 
(Clark et al. 2008); it has an equatorially concentrated mass 
of hot (T ^ 1000 K) dnst having a radius of 60-600 AU; 
and it exhibits a similar correlation between emission-line 
width and ionization potential. The presence of strong He ll 
emission in NaStl, however, suggests a significantly hotter 
ionization source than that of Wdl-9, which does not exhibit 
Hell emission. 

Clark et al. (2013) interpreted the double-peaked line 
profiles of Wdl-9 as being the result of quasi-Keplerian ro¬ 
tation in a massive disk or torus. This interpretation has 
also been invoked to explain the double-peaked line profiles 
detected in the LBV binary MWC 314 (Lobel et al. 2013; 
also shown in Figure[9]). A rotational component from the 
inner regions of NaStTs disk is thus worth considering as 
a potential explanation for the double-peaked profile. How¬ 
ever, such profiles can also be produced by another mecha¬ 
nism. An expanding disk, torus, or flattened hollow cone ge¬ 
ometry can generate double-peaked line profiles, as demon¬ 
strated by Zickgraf (2003), who argned that many sgB[e] sys¬ 
tems could actually be exhibiting a combination of rotating 
disk emission at relatively small radii (where permitted lines 
are formed) and emission from a radially outflowing disk at 
larger radii (where forbidden lines, like [Nll], are formed). It 
is difficult to discriminate between these possibilities for the 
double-peaked profiles from the inner regions of NaStl, al¬ 
though the hypothesis of quasi-Keplerian rotation might be 
more compatible with the observed correlation between pro¬ 
file width and ionization potential shown by CS99. On the 
other hand, the kinematics of the outer nebula, as indicated 
by Figure[4l are more likely to be the result of an equatorial 
outflow that is not perfectly azimuthally symmetric, since 
rotation at such large radii seems unlikely. 

5.2.1 NaStl as an advanced analog to RY Scuti 

Another very relevant object for comparison to NaStl is 
the eclipsing binary RY Scuti. It consists of a massive 
O9/B0 supergiant primary currently undergoing RLOF, and 
an accreting secondary (likely to be an 05 star) that is 
buried inside an opaque accretion torus (Grundstrom et al. 
2007; Smith et al. 2011b). The system as a whole is sur¬ 
rounded by an expanding double-ring nebula of equatori¬ 
ally concentrated ionized gas exhibiting multi-peaked emis¬ 
sion profiles from both Ha and [Nil] (Smith et al. 2002), 
the latter of which is also shown in Figure[9l In this case, 
the multi-peaked morphology is clearly the result of non- 
azimuthally symmetric expansion, not rotation. The rings 
are «1800 AU in extent, expanding at roughly ±45kms“^. 
Based on proper-motion measurements, the structure ap¬ 
pears to have formed from multiple discrete eruptions of 
mass that occurred 130 and 250 yr ago. Smith et al. (2011b) 
suggested the possibility that the equatorial outflow was pro¬ 
duced by an LBV-like stellar eruption that occurred in re¬ 
sponse to the sudden accretion of mass, and that the sys¬ 
tem is potentially on its way to becoming a B[e] supergiant 
-I- stripped-envelope WR binary. Interpreting NaStl in light 
of RY Scuti is interesting because NaStl already exhibits 
evidence consistent with the presence of both a B[e]-like 
or 0[e]-like supergiant star (producing the double-peaked 


emission-line profiles) and a hot WR (to provide the ioniza¬ 
tion and, perhaps, the X-ray emission), all surrounded by 
an equatorially-concentrated nebula of CNO-processed ma¬ 
terial having a physical extent over an order of magnitude 
larger than the ionized nebula around RY Scuti. It is thus 
plausible that NaStl is a more highly evolved analog of the 
RY Scuti system. 

5.3 On the lack of WR spectral features in NaStl 

CS99 concluded that the central ionization source in NaStl 
has a stellar temperature in excess of 30 kK, and is thus 
likely to be an early-type WR. Such objects typically ex¬ 
hibit fast winds that generate strong emission lines with 
FWHM velocities of « 2000-3000 kms“^. The absence of 
any such broad spectral lines from the purported early-type 
WR is thus puzzling; we speculate below on several potential 
causes. 

Heavy obscuration of the WR spectrum by dust could 
be exacerbated by the inclination of the system with respect 
to our line of sight (« 12° from edge on). Dust formed via 
colliding winds or in an equatorial outflow will probably be 
concentrated in a similarly inclined plane, which, if suffi¬ 
ciently compact, could obscure the central stars, while still 
allowing some photoionizing radiation to escape at higher 
latitudes. We speculated in §3 that the resolved structure in 
the Ks band could be the result of Hel A2.06/rm emission 
from the central source or stellar wind scattered beyond the 
obscuring nebula, which would be allowable for the case of 
an equatorial/toroidal distribution of dust. However, if dust 
obscuration is to blame for the heavy suppression of the stel¬ 
lar spectral features, then it is difficult to explain why the 
same dust does not also obscure the high-ionization double- 
peaked emission lines detected by CS99 or the 1-2 keV X-ray 
emission shown in §4, which should arise in the hot inner re¬ 
gions of the system, very close to the central stars. 

An alternative explanation for the lack of WR fea¬ 
tures is dilution by a more luminous companion star and 
the circumstellar material that produces the bright nebular 
spectrum. After all, WRs that form via mass-transfer can 
stem from stars that have significantly lower initial masses 
(> 8 Mq; Podsiadlowski, Joss, & Tsu 1992; Vanbeveren et 
al. 1998) and lower luminosities than WRs that form from 
direct mass loss via strong winds (> 20-25 M©; Meynet 
& Maeder 2005; Ekstrom et al. 2012; Smith & Tomble- 
son 2015). Lower-mass WRs will likely have weaker winds 
(slower and less dense) than classical WRs of relatively high 
mass and luminosity. As such, these stars might not be 
easily recognizable spectroscopically as WRs, lacking sim¬ 
ilarly broad and luminous emission lines that can easily 
shine through dense circumstellar material, as in the case 
of NaStl’s nebula. Relatively slow winds were already noted 
as a possible explanation for the relative faintness of NaStl 
at harder X-ray energies relative to many other WR CWBs. 
It thus seems plausible that dilution of the WR spectrum 
by a luminous companion and nebula, combined with some 
obscuration of the WR by dust, is the reason for the lack of 
detectable WR features in the optical spectrum. Similarly, 
the infrared spectrum of the WR could be diluted by the ex¬ 
tremely luminous (it's ~ 6.3 mag) hot dust component. The 
K-band brightness of a relatively under-luminous, reddened 
WR at the distance of NaStl (adopting Mk = —3 mag for 
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faint WRs from Crowther et al. 2006) is expected to be ~ 2- 
4 mag fainter than the thermal dust, even without account¬ 
ing for addition circumstellar extinction, so heavy dilution 
of the stellar spectrum in the infrared is also plausible. 


6 SUMMARY AND CONCLUDING REMARKS 

Based on the results of high-resolution optical imaging with 
HST, X-ray imaging with Chandra, and AO infrared imag¬ 
ing with Magellan, and the nebular analysis of CS99, we 
conclude that NaStl (WR122) is likely to be a WR binary 
enshrouded in an equatorial nebula of CNO-processed gas 
with a compact thermal dust component. The apparently 
slow kinematics and equatorial geometry of the nebula sug¬ 
gest that it is the byproduct of mass transfer. NaStl could 
thus be an analog to massive interacting binary systems, 
such as RY Scuti, MWC 314, and Wdl-9, but at a more 
advanced evolutionary stage, where one of the components 
has lost enough mass to become a WR. NaStl could there¬ 
fore represent a rare and important example of the recent 
formation of a stripped-envelope WR formed from binary 
interaction. 

Many open questions about the NaStl system remain, 
particularly regarding the origin of various spectral features, 
dust chemistry, the physical parameters of the component 
stars, and the orbital period of the system. Such information 
might be obtainable via long-term photometric and monitor¬ 
ing and high-resolution spectroscopy in the optical and in¬ 
frared. Detailed spatially-resolved chemical abundances will 
also shed more light on the recent evolution of the compo¬ 
nent stars and the processes that created the nebula. 
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